Magnetic Order and Spin Dynamics in the Proximity of a Ferromagnetic Quantum 

Critical Point: a fiSR study of YbNi 4 P 2 
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The local 4/-electronic spin dynamics and magnetic order in YDM4P2 were studied by means of 
muon-spin relaxation measurements. Zero-field muon-spin relaxation proves static magnetic order 
with a strongly reduced ordered Yb 3+ moment of (2.5 — 4.6) x 10~ 2 /ib, below Tc = 140 mK. 
Above Tc, the muon spin polarization P(t, B) is dominated by quasihomogeneous spin fluctuations 
and exhibits a time-field scaling relation P(t,B) = P(i/£> 7 ), indicating cooperative critical spin 
dynamics in the system. At T = 190 mK, slightly above Tc, 7 = 0.81(5) suggesting time-scale 
invariant power-law behavior for the dynamic electronic spin-spin autocorrelation function. 

PACS numbers: 71.27.+a, 75.30.-m, 76.75.+i, 75.50.Cc 
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Lanthanide-based heavy-fermion (HF) systems are 
suitable model systems to study emergent phenomena at 
a quantum critical point (QCP), where collective quan- 
tum fluctuations trigger the system continuously from a 
magnetically ordered to a non-magnetic ground state 1 - - — . 
However, despite intense research, to the best of our 
knowledge, no 4/-based material is known with a contin- 
uous ferromagnetic (FM) to paramagnetic (PM) quan- 
tum phase transition (QPT). The existence of such a 
QPT is also controversially discussed from a theoretical 
point of viewer—. 

Recently, Krellner et al suggested that the HF metal 
YbNi4P2 with a quasi-one-dimcnsional (1-D) electronic 
structure exhibits FM quantum criticality above a low 
FM transition temperature T c = 170 mK 12 . YbNi 4 P 2 
crystallizes in the tetragonal ZrFe4Si2 structure contain- 
ing isolated chains of edge-connected Ni tetrahedra along 
the c— axis. The Yb atoms are located in the channels 
between these Ni tetrahedral chains. The reduced dimen- 
sionality in the Yb and Ni network and the geometrical 
frustration between neighboring Yb chains give rise to en- 
hanced quantum spin fluctuations of the magnetic Yb 3+ 
ions. In the PM state above 50 K, the magnetic sus- 
ceptibility shows Curie- Weiss behavior with an effective 
moment fJ- e ff = 4.52/is that is characteristic for mag- 
netic Yb 3+ ions. Analysis of the magnetic entropy re- 
veals a Kondo energy scale of Tk ~ 8 K for the crys- 
tal electric field ground state doublet. The FM transi- 
tion is evidenced by distinct anomalies in magnetic sus- 
ceptibility, specific heat, and resistivity measurements. 
Low-T magnetization measurements suggest an ordered 
FM moment of m or d ~ 0.05(4)/ig. Pronounced non- 
Fermi-liquid (NFL) behavior is reflected by a stronger- 
than-logarithmic diverging Sommerfeld coefficient and a 
linear-in-T resistivity state apparent in a T range larger 
than a decade above Tc- In external magnetic fields, 
the NFL behavior is suppressed and FL behavior grad- 
ually recovers. Therefore, YbNi4P2 is considered as a 



clean system situated in the very close vicinity of a FM 
QCP, with FM quantum fluctuations dominating ther- 
modynamic and transport quantities at T > Tc- 

The present knowledge on YbNi4P2 is based on mea- 
surements of macroscopic magnetic, thermodynamic, and 
transport properties. The next step in a deeper investi- 
gation of this prospective FM quantum critical system 
is to get insight on a microscopic level. Beside the na- 
ture of the magnetic order, a central issue in the present 
context of critical behavior is the spin dynamics. Since 
in systems close to a QCP, the ordered moment is usu- 
ally strongly reduced, muon spin relaxation (/iSR) has 
proven to be an extremely valuable technique to collect 
appropriate information 1 ^—. 

Here, we present /iSR experiments on polycrystalline 
YbNi4P2, providing microscopic evidence for static mag- 
netism at T < Tc ~ 140 mK with an ordered mo- 
ment of m or d — (2.5 — 4.6) x 10~ 2 ^s/Yb, depending 
on the assumed muon site. Above Tc, the muon-spin 
polarization P(t) obeys the time-field scaling relation 
P(t) = P(t/B - 81 ^), indicating cooperative and critical 
spin dynamics. 

In a /zSR experiment positive spin-polarized muons are 
implanted into the sample, and the subsequent time evo- 
lution of the muon spin polarization is monitored by de- 
tecting the asymmetric spatial distribution of positrons 
emitted from the muon decayi^. /iSR in longitudinal ap- 
plied magnetic fields is dominated by Yb-4/ electronic 
spin fluctuations that couple to the implanted muons. 
The /LtSR experiments on YbNi 4 P2 in zero field (ZF) 
and longitudinal (LF) applied field - with respect to the 
initial muon spin polarization - were performed on the 
7rM3 beam line at the Swiss Muon Source (S/iS) at the 
Paul-Scherrer-Institut, Switzerland. The sample was pre- 
pared by crushing ~ 270 mg of single crystalline mate- 
rial, grown in a self-flux at 1400 °C in a closed Tantal 
crucible and characterized by powder x-ray diffraction ex- 
periments, proving the absence of any foreign phases. De- 
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FIG. 1: (Color online) (a) Corrected muon spin polariza- 
tion P(t) at ZF for representative T above and below Tc ~ 
140 mK. At T < Tc, solid lines are fitting curves accord- 
ing to Eq. (1). (b) T dependence of the ZF fiSH frequency 
/ M (T). The solid line is a fit to the phenomenological function: 
U = Ui°) ■ C 1 - W^T- ( C ) T dependence of the ZF static in- 
ternal field distribution a in Eq. (1). The solid line is a guide 
to the eye. (d) T-dependence of 1/TiT in the PM regime. 

The line describes power-law behavior as — oc T~ . 
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tailed low-T measurements on polycrystalline YbNi4P2 
were reported elsewhere^. 

Figure 1(a) displays typical time dependencies of the 
ZF muon-spin polarization P(t) in YbNi4P2 at repre- 
sentative temperatures above and below Tc- A finite 
T-independent background signal due to muons that 
stopped in a Ag sample holder (signal fraction w 50%) 
was taken into account. At T > 160 mK, an exponen- 
tial muon-spin relaxation is associated with fast fluctu- 
ating paramagnetic electron spins with a relaxation rate 
A(160 mK)« 0.152(2)/is _1 . Note, that dense static nu- 
clear dipole moments would give rise to a weak Gaussian 
relaxation in the PM regime. While cooling through Tc, 
an additional magnetic relaxation mechanism is appar- 
ent, strongly increasing with lowering T. Below Tc, & 
low-frequency oscillation with a Gaussian relaxation of 
the muon-spin polarization is observed indicating mag- 
netic ordering of weak electronic Yb 3+ moments. The 
muon-spin asymmetry data in the FM regime can be de- 
scribed best using the functional form 17 ' 18 - 



Pit) = \ + ^[cos(27r/^) - • sin(2^t)] • e^' 2 ,(1) 

where /„ and a are the muon spin precession fre- 
quency and the Gaussian field width, respectively. The 
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FIG. 2: (Color online) Corrected muon spin polarization at 
T — 20 mK and various longitudinal magnetic fields Blf- 
The lines represent theoretical depolarization curves for the 
static GKT function in corresponding longitudinal fields. 



2/3 oscillating and the 1/3 non-oscillating terms origi- 
nate from the spatial averaging in polycrystalline sam- 
ples, where 2/3 (1/3) of the internal magnetic field com- 
ponents are directed perpendicular (parallel) to the ini- 
tial muon spin, causing a precession (no precession) of 
the muon spin. The observation of a 2/3 and 1/3 sig- 
nal fraction below Tc, implies dense magnetic moments 
and proof that 100% of the sample volume shows static 
magnetic order. The latter is supported by LF-/iSR mea- 
surements as discussed in detail below. In the limit 
2tt//j 3> cr, Eq. (1) becomes a Gaussian damped cosine 
function. For 27r/ M — > 0, close to the magnetic transi- 
tion, Eq. (1) is equivalent to the Gaussian Kubo-Toyabe 
(GKT) function^, which describes a muon-spin relax- 
ation due to a static Gaussian field distribution centered 
around Bi oca i — 0. ZF-//SR on the antiferromagnetically 
ordered system YbRh^Si^ reveals a similar crossover 
from a Lorentzian to a Gaussian damped /iSR signal 
in the vicinity of the PM to magnetic phase transition, 
attributed to a transition from dynamic to static mag- 
netism of magnetic Yb 3+ moments. 

For YbNi4P2 a finite fj.SK frequency is clearly observed 
below 150 mK. From the measured frequency value / M = 
0.188(1) MHz at 20 mK one can determine the internal 
local field at the muon site to Bi oca i = 13.87 G using 
Biocai = 271-/^/7^ with 7 M = 2tt x 13.55 kHz/G as the 
muon gyromagnetic ratio. The local field Bi oca i as well 
as the local static field width ABi oca i = cr/7 M , » 6 G are 
very small for conventional rare-earth magnets with large 
ordered moments. The fractional width ABi oca i/ Bi oca i 
of the spontaneous field distribution is ~ 0.4 at low T 
and remains constant as T -> Tc, which is a reason- 
able value for a magnetically ordered HF system, as e.g. 
in CeRhIn 5 ABi oca i/Bi oca i = 0.5 is observed^. Thus, 
the local field distribution is nearly uniform and homo- 
geneous in the FM regime. The spontaneous muon-spin 
precession and the Gaussian shape of the internal field 
distribution below Tc arise from a dense system of weak 
magnetic moments with small, static magnetic inhomo- 
geneities. The presence of a finite Bi oca i ^ proves co- 
herent magnetic order. 
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FIG. 3: (Color online) Main panel: Field dependence of the 
dynamic muon spin relaxation rate Al- The solid curve rep- 
resents a Redfield fit. For display reasons, the ZF value is set 
at Blf = 0.01 G. Inset: Field dependence of the corrected 
muon spin polarization P(t) at T = 190 mK. 
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FIG. 4: (Color online) Corrected muon decay asymmetry at 
T = 190 mK for various magnetic fields as function of the 
scaling variable t/B^'p 1 . The dashed-dotted line is a fit of the 
13 G data with P(t)/P(0) = exp[-A L t]" a9 . 



ZF-/iSR allows a precise determination of the T de- 
pendence of the magnetic order parameter, which is pro- 
portional to the measured /xSR frequency / M . The T 
dependence of / M and a is shown in Figs. 1(b) and 1(c). 
For T < 140 mK, both observables exhibit a continu- 
ous increase. The T dependence of / M can be fit to the 
phenomenological function = / M (0) ■ (1 — ^)™ f° r 
T < T c with n = 0.208 ± 0.02, / M (T) = 0.199(3^ MHz, 
and Tc = 140(2) mK. The value of the effective critical 
exponent n, describing the critical behavior close to Tc, 
is between n = 0.125 and 0.325, which are theoretically 
expected for two-dimensional (2D) and isotropic three- 
dimensional (3D) Ising magnets, respectively. This is not 
in contradiction with the claim of a quasi- ID system. In 
such a system, the weak inter-chain coupling results in 
an evolution from a ID behavior at high T to a 2D Ising 
or 3D behavior at low T, which is intimately linked with 
(and is a prerequisite for) the long-range ordering at finite 
T. The low data point density between 0.6 < ^ < 1, 
however, precludes the determination of the precise criti- 
cal exponent. The obtained value for Tc agrees well with 
the value found in specific heat measurements on these 
single crystals^. 

For all examined T < Tc, the sample signal is an- 
alyzed with a well-defined single / M and a, signaling 
that the magnetic order is a bulk effect and that only 
one dominant muon stopping site is present. In gen- 
eral, for the determination of the muon stopping site(s) 
it is important to deduce the hyperfine coupling con- 
stant. One way to find potential muon sites is to com- 
pare calculated and measured quantities for the local field 
Biocai at the muon site. The muon preferentially settles 
at tetrahedra or octahedra interstitial crystallographic 
sites. From simple symmetry arguments the most prob- 
able muon stopping sites, using WyckofFs notation, 
are 4/(1/4,1/4,0), 8j(l/4,l/4,l/4), 4/(1/4,1/4,1/2), 
8i(l/4,l/2,l/2), 4c(l/2,0,0), 4c(l/2,0,l/2), 26(0,0,1/2), 
and 2a(l/2,l/2,l/2). For a particular FM structure with 
the magnetic Yb 3+ moments aligned within the a b plane 



and a dominant 4/-/i dipolar interaction, one can de- 
termine the expected internal field values for the pro- 
posed sites. Our lattice sum calculations reveal that only 
at the 4c(l/2,0,l/2) site and 8i(l/4,l/2,l/2) site a lo- 
cal field Biocai of the measured absolute magnitude is 
found. For the 4c and 8i sites the measured local field of 
Biocai — 13.87 G corresponds to a static ordered moment 
of the Yb ions of m or d = 0.046/1 s and 0.025/is, respec- 
tively. Both values are in good agreement with the value 
deduced from recent magnetization measurements 12 . 

The temperature dependence of the exponential relax- 
ation rate A/, = y-, observed above Tc, is plotted in 
Fig. 1(d) on a log- log scale as Cooling down from 

800 mK, exhibits power-law behavior according to 
_i_ x r -i.40(6)_ At T < 190 mK, the power-law behav- 
ior in jhp persists in the PM regime down to Tc, how- 
ever, with a slight change of the critical exponent, i.e., 
7pkp oc T -1 - 5 W (dashed-dotted line). The observed 

behavior is close to the T~ 4 / 3 temperature dependence 
predicted by the self-consistent renormalization (SCR) 
theory for a system close to a 3D ferromagnetic QCP 22 . 
There is no prediction for an itinerant quasi ID system in 
the T range between the exchange energy scale and order- 
ing temperature. For an isolating ferromagnetic quasi ID 
spin chain the T dependence of the relaxation rate above 
Tc depends strongly on the details of the interactions - 
see, e.g.r^. 

LF-/iSR experiments allow to separate the dynamic 
contribution to the relaxation of the muon-spin polar- 
ization. Investigations of the low-T muon-spin dynamics 
yield additional information about the origin of the NFL 
behavior in YbNi4P2. Figure 2 displays the muon-spin 
asymmetry function P(t) at T — 20 mK for different 
applied LF's. The muon-spin relaxation is completely 
suppressed in an applied field Blf ~ 300 G, demonstrat- 
ing that the internal field distribution is static in nature. 
However, the observed decoupling can not be described 
accurately by a standard muon asymmetry function that 
considers an internal field distribution which is symmet- 
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ric around Bi oca i = 0. For comparison, Fig 2. shows the- 
oretical depolarization curves for the static GKT function 
in the corresponding longitudinal magnetic fields. This 
supports the ZF data, i.e., the observation of a broad field 
distribution centered around a finite but small internal 
field S /oca/ (20 mK) w 13.87 G in the FM phase. Finally, 
when Blf 3> Bi oca i, the muon spin relaxation is decou- 
pled from the static Bi oca i as observed for Blf > 23 G. 

At T > Tc, the field dependence of the muon-spin 
relaxation probes the Fourier transform of the dynamic 
spin-spin autocorrelation function q(t) — (Sj(t) • Sj(0)), 
which exhibits exponential behavior for homogeneous 
systems and power-law (or cutoff power-law) or stretched 
exponential behavior for inhomogeneous systems. The 
inset of Fig. 3 displays the muon-spin polarization P(t) 
at T = 190 mK, both in magnetic LF between 13 and 
143 G and ZF. The relaxation rate Xl is reduced with 
increasing field. The field dependence of Xl is given in 
the main panel of Fig. 3. It shows nearly no field de- 
pendence for magnetic fields of less than ~ 13 G, but 
varies more strongly, as H~ K with k ~ 0.79(7), for higher 
fields. From the field dependence of Xl, the spin autocor- 
relation time t c can be estimated using the Redfield for- 
malisrn for X L (B LF ) = (2^(B} 1uc )t c )/[1 + ( 7 2 5| f t c 2 )] 
considering r c as independent of the applied field Blf- 
Here, Bfi uc (t) describes the time-varying local magnetic 
field at the muon site due to fluctuations of neighboring 
Yb 3+ moments, with a local time averaged second mo- 
ment A 2 = 7 2 (-B|; uc ) and a single fluctuation time r c . 
For Tvjj <C fcflT (ui giving the spin fluctuation rate), the 
fluctuation-dissipation theorem^! relates r c to the imagi- 
nary component of the local g-independent /-electron dy- 
namic susceptibility, i.e. T C (B) = (kBT)\x"(w)/u)]. The 
fit to the data (solid curve in the main panel of Fig. 3) 
yields A 2 » 0.1 (MHz) and r c sa 6 x 10~ 7 s, the latter 
value nearly three orders of magnitude larger than the 
one obtained for YbRli2Si2 at T — 20 mK' 4 :, suggesting 
very slow critical fluctuations. 

The /iSR time spectra in Fig. 3 are well described with 
a stretched exponential relaxation function of the form 
P(t) = P(0) exp[-(Ai)' 3 ]. An exponent of « 0.9 shows 
that the relaxation rate is nearly uniform throughout the 
sample, indicating that YbNi4P2 exhibits quasihomoge- 
neous spin fluctuations for T <C Tr. The spin dynam- 
ics is characterized by a narrow distribution of correla- 
tion times {(3 — 1 corresponds to one single correlation 
time). Thus, disorder-driven theories, including Kondo 
disorde r 25 ' 26 and the Griffith phase scenario 27 as primary 
mechanisms for the observed NFL behavior, can be ruled 
out. It further implies that the crystalline disorder in 
YbNi4P2 is quite small, which is consistent with a small 
residual resistivity (po ~ 2.4/if2 cm) and the stoichio- 
metric occupation of the crystallographic lattice sites re- 
vealed by the x-ray structure refinement^. 

A sensitive test to identify power-law or stretched ex- 



ponential behavior of q(t) is a time- field scaling analysis 
of the muon-spin relaxation function. In both cases a 
specific time-field scaling can be found, i.e., the muon- 
spin relaxation function P(t,BLF) obeys the scaling re- 
lation P(t,BLF) = P(t/Bl F ). This relation applies 
only in the asymptotic strong field limit, i.e., as long as 
27r/ M = j^Blf 3> Xl 28 - If time-field scaling is obeyed, 
a plot of P(t, Blf) versus t/B^ F at T > T c will be 
universal for the correct choice of 7, and distinguishes 
between power-law (7 < 1) and stretched exponential 
(7 > 1) correlations. For small Blf, the field depen- 
dence is expected to be due to the change of / M rather 
than an effect of field on q(t). A breakdown of time-field 
scaling would occur for high fields where q(t) is directly 
effected by the applied fields. Figure 4 shows the same 
asymmetry data, as displayed in Fig. 3, as functions of 
the scaling variable t/B'[ F . For 7 = 0.81(5) the data 
scale well over ~ 2.5 orders of magnitude in t/ B LF and 
for all fields between 13 and 143 G, except for 293 G. 
Here, at large t, the data fall above the low-field scaling 
curve. Fields /ib-Blf > ksT (with ks =Boltzmann's 
constant) would be expected to affect the spin dynam- 
ics. The scaling exponent 7 = 0.81(5) < 1 implies that 
within the /iSR frequency range, the spin-spin correla- 
tion function q(t) is approximated by a power law (or a 
cutoff power law) rather than a stretched exponential or 
exponential 28 , consistent with the Redfield analysis. The 
power-law is time-scale invariant and dynamical modu- 
lations should therefore be observable in any time win- 
dow. The obtained time-field scaling of the relaxation 
data is a signature of slow homogeneous spin dynamics. 
It strongly indicates that the critical slowing down of 
spin fluctuations at the magnetic phase transition occurs 
cooperatively throughout the sample. In stoichiometric, 
homogeneous NFL systems such behavior may arise from 
the effect of disorder on quantum critical fluctuations in- 
herent to a QCP. This is suggested for the NFL com- 
pound YbRhaSi?^ 5 -. 

In conclusion, ZF-/iSR in the stoichiometric NFL com- 
pound YbNi4P2 clearly proves static magnetic ordering 
of strongly reduced Yb 3+ moments below Tc — 140 mK. 
Above Tc, the muon spin polarization P(t) obeys the 
time-field scaling relation P(t) = P(t/5°- 81 ( 5 )) for ap- 
plied magnetic fields B between 13 and 143 G, indicat- 
ing cooperative and critical spin dynamics. Power-law 
behavior of the dynamic spin-spin autocorrelation func- 
tion is implied by the observation of 7 < l 28 . The LF- 
/iSR results suggest that the NFL behavior observed at 
T > Tc is induced by quasi homogeneous critical spin 
fluctuations. 
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